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£S) ■ Abstract 

In this article we apply a recently invented analytical real-space renormaliza- 
tion group formulation which is based on numerical concepts of the density ma- 
fS^ trix renormalization group. Within a rigorous mathematical framework we con- 

struct non-perturbative renormalization group transformations for the spin-i XXX 
Heisenberg model in the finite temperature regime. The developed renormalization 
group scheme allows for calculating the renormalization group flow behaviour in the 
temperature dependent coupling constant. The constructed renormalization group 
transformations are applied within the ferromagnetic and the anti-ferromagnetic 
regime of the Heisenberg chain. The ferromagnetic fixed point is computed and 
compared to results derived by other techniques. 
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I. INTRODUCTION 



In 1966 L.P. Kadanoflu presented arguments which would allow one to calculate critical ex= 
ponents without ever working out the partition function explicitly. Later in 1971 K.G. Wilsona 
transformed Kadanoff's qualitative block spin explanations into a quantitative formulation in the 
context of critical phenomena which has become known as the Wilson Renormalization Group 
(RG) technique. 

The core physical concept of the RG is the scale invariance at the critical point. As the phys- 
ical system moves towards a phase transition, it becomes increasingly dominated by large-scale 
fluctuations. At the critical point the correlation length, i.e. the length scale of the fluctuations, 
becomes infinite and the system exhibits scale invariance. Within an application of the RG a 
RG transformation (RGT) needs to be defined which eliminates non-relevant degrees of freedom 
of the system. If the system becomes invariant by carrying out successive RG transformations, 
a RG fixed point is reached. A statistical system exhibits two trivial fixed points, belonging to 
zero and infinite temperature, where the system possesses inherent scale invariance according to 
complete order or complete disorder respectively. A nontrivial fixed point, if present, corresponds 
to a critical point of the physical system. The behaviour of physical quantities at the critical point 
of the system is described by scaling laws using critical exponents. 

Although RG methods have been successfully applied to a variety of physical problems, the con- 
struction of RGTs applicable to strong coupling regimes is, apart from a few exceptions, an un- 
K> , solved and challenging problem. ExampkjjLof current research are the strong coupling quantum sain 

chains, including the Heisenberg modelsLro, and nonlinear partial differential equations (PDEs)ErQ. 
Standard approaches like perturbation theory in combination with Fourier space RG methods can- 
not be applied. 

In this paper we explore the critical behaviour of the isotropic spin-i Heisenberg model at finite 
temperature using the Generalized Real-space Renormalization Group (GRRG). The method was 
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introduced in an earlier work as the general (real-space) RGLJ. The GRRG requires the definition 
of an auxiliary space 7i a ux prior to the construction of the desired local GRRG transformation 
(GRRGT), i.e. the RGT for a subsystem called a block, analogous to the Kadanoff block spin 
formulation. The auxiliary space describes the quantum correlations resulting from the different 
boundary conditions between separated adjacent blocks within the local GRRGT. If the defined 
auxiliary space allows for a decomposition of the quantum system into blocks keeping all possi- 
ble boundary conditions the local RGT is called exact. In this work we construct a perfect local 
GRRGT based on an auxiliary space providing an approximate description of the boundary con- 
ditions. The local GRRGT is formulated as a composition of two linear maps, the embedding 
map G and the truncation map G + . Both maps depend on the choice of the auxiliary space and 
are constructed according to an imposed physical constraint, the invariance relation!]. Originally 
the conser|iation of the free energy or the partition function was used to formulate the physical 
constraint .0. _ 

In this work we use the notation of the original worku. Block quantities are indexed by capital 
letters, corresponding to sites in the blocked chain. The indexing set for the blocks is denoted 
as 3. Neighbouring blocks are indexed by a sequence I, I — 1, J — 2, • • • € 3 whereas independent 
blocks are indexed by different letters I, J, ■ ■ ■ £ 3. A block Hilbert space Hi contains at mini- 
mum two single site Hilbert spaces TLi and TU-i- Single site Hilbert spaces are indexed by letters 
i, j,k, . . . . If a single site space TLi is contained within a block Hilbert space TLi we write TLi C TLi 
or i € I if it is obvious that / refers to the block Hilbert space. We further use the abbreviation 
{i,i — 1,. ..} C {1,1 — 1,. ..} instead of W{j ) j_i ! ...} C TL{iy (g> TL{i-i}, .... Using this notation it 
is not apparent which single site space is contained in a particular block Hilbert space. If this is 
important it needs to be pointed out explicitly. 

In the next section we begin by revisiting the classical case. No quantum correlations occur in the 
classical analogue of the Heisenberg chain and the concept of an auxiliary space is therefore un- 
necessary. However, we construct a local RGT and by direct comparison with the GRRG method 
we provide the reader with a feel for the abstract mathematical formulation of the GRRG method. 



We proceed in section [II by constructing a perfect local RGT for the spin-^ XXX Heisenberg 



model following the concepts in the original workQ. In section IV we discuss analytical results 
of the flow behaviour calculated by the perfect GRRGT. Critical exponents are calculated from 
the nontrivial ferromagnetic RG flow behaviour for the three dimensional Heisenberg chain and 
compared with results calculated by other methods. In section [V] we derive the flow behaviour for 
the anti-ferromagnetic regime of the one-dimensional Heisenberg model. In the final section we 
conclude with some perspectives on the GRRG method. 



II. THE MIGDAL-KADANOFF RGT 

We consider the one dimensional Ising model without an external magnetic field and with nearest 
neighbour (n.n.) interactiontl. All thermodynamic quantities can be calculated from the model's 
partition function 

N 

Z Ising =Y, e ~ 1/(kBT)HlBias with ff Ising = J2 Ja i0 - i+1 , (1) 

where {<7i} denotes that the sum should be extended over all possible assignments of ±1 to each 
lattice site i corresponding to an array of elementary spins {(Ji\ placed on the lattice sites {i}. In 
(pi) ks is Boltzmann's constant and we are interested in the limit N — ► oo. 
Introducing a temperature dependent coupling constant 

K = ST (2) 



we write the partition function in the form 

N/ 2 
-Rising = X! IT e ^Pi K (^i^l^i + <J2iCr2i+l)} ■ (3) 

Decomposing the sum over all possible configurations into odd and even sites the sum over the 
even sites is calculated by successive application of 

^ exp [K ((T2i-i<J2i + o-2jCT2J+i)] = exp [K (er 2 j-i + v%i+i)\ + exp [-K (a 2 i-i + °2i+i)] (4) 

0-2* =±i 

for every even site 2i. The application of the GRRG method requires a physical constraint to 
construct a local GRRGT imposed by an invariance relatione. We define the classical analogue of 
the invariance relation by keeping the partition function unchanged 

Z[H lsiag (J)} = Z[H lsing (J'J')} . (5) 

In equation (H) /' denotes a change in the ground-state energy of the energy function Hi s j ng (J). 
By inserting (pi) and (0) into equation (ph we define the effective functional dependence as 

N/ 4 

0{{j,}, K', /') = J| exp [K 1 {a 2l ^a 2l + a 2l a 2l+1 ) + 4f/N] . (6) 

i 

Using the definition (|6|) we derive the classical analogue of the global GRRGT as 

N/2 N/4 

^ Yl exp [K ((J2i-\<Tii + O2»02;+i)] = ]~J exp[K' (a 2l -icr 2i + <T 2i <T 2i+ i) + Af'/N] . 

{a-j , j even} i i 

(7) 

The GRRG method requires a decomposition of the spin chain into commuting blocks, which can 
always be performed in the classical case. We therefore write the local GRRGT as 

^ exp[K(a 2l ^ 1 (j 2i + <J 2 icr 2 i+i)} = exp [K 1 a 2l -i <r 2i+1 + if/N] (8) 

CT 2i =±l 

and the effective coupling is calculated as K' = (1/2) • lncosh(2iir) which yields-tp the trivial RG 
flow behaviour as it is expected for one dimensional strongly correlated systema!3. Relation (||) is 
the classical analogue of an exact local GRRGT since the invariance relation (g) can be derived 
from the local GRRGT. 

The previous calculations are a reinterpretation of the Migdal-K adan off transformation for classical 
spin systems. The calculation was fixst done by A. A. MigdalEj-O and reformulated using bond 
moving techniques by L.P. KadanoflEm. 



III. THE CONSTRUCTION OF THE LOCAL RGT FOR THE ISOTROPIC 

HEISENBERG CHAIN 

In this section weJerive a perfect local RGT for the isotropic spin-^ Heisenberg chain by applying 
the GRRG methodQ. To the best knowledge of the author no other controllable approximation is 
currently available to analytically calculate the critical properties for the quantum spin chain at 



finite temperature. 

The Hamiltonian of the model is defined by 



H 



XXX 



N 



(9) 



which is totally isotropic in the spin components and known as the XXX spin-i modeltSlia. The 
spin variables a x , a y and a z define the Lie algebra st(2). In this article we choose the smallest 
nontrivial representation Sf = (H/2)af by the Pauli matrices 
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The partition function for the one-dimensional Heisenbcrg model is defined by 



(10) 



-Zxxx =tr {CT3 .} exp 



te{a 3 } exp 
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(11) 



where we used the vector notation and introduced a temperature dependent coupling constant 
K = ~ T - In equation ( |l l| ) trr CT .i denotes the trace over all lattice sites in the quantum chain. 
Analogous to the classical case reported in section O the partition function is used for defining the 
invariance relation 



-Zxxx [O ({a j },K)] = tr {crj 1 3 odd} tr {CTj _ i cvon} exp 



N 



Kj~](Ti- a i+1 



i=l 



tr 



{o-j, j odd} ex P 



N/2 



f'{K) + K\K) J2 *m-i • ffw+i 



j=i 



Z xxx [0({a 3 , j odd}, K'J')} 



(12) 



where we used of the factorization property of the trace tr r^.i = J^ t r o-;- 

According to the construction of the local GRRGT -ase proceed by changing our notation and 
equip every operator with an abstract auxiliary spaced which is currently not further specified. 
The action of each of the operators on the auxiliary space is defined as the identity map until 
further specifications are given. The embedding map Gn'®H'^ x an d the truncation map G^ 8W 
are defined according to the invariance relation 

-Zxxx [Oh«h„ x (K)] = Zxxx \pH®n ° ^§h„ x (K) o GH®H aux ] = Zxxx [Owgw^ (K')j • 

(13) 

The embedding map Gu'mH 1 together with the truncation map G^^ w define the GRRGT. 
By choosing 
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H®H 



aux (K) = exp 
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k y^ Si ■ (Tj+i 



with K = (K, f) 



(14) 



the embedding and truncation maps are defined as 

G W8Mmx = tr{ ffj; jovcn} ® l{ ( 7 i ,jodd} <8 l« aux and Gh'®h'.^ = Ih'qH'^ ■ (15) 

Analogous to the classical case the local operators are derived from their global counterparts by 
decomposing the quantum chain into blocks. Taking the trace over one even site in each block we 
arrive at the relation 



N 
-ZxXX [CWiWaux (K)] = tr^joddjtr^jeven} exp 



■K" / , <?i • <?,:+i 



i=l 



® 1 



Ha 



N/2 



tr^.jodd} JJ trcvonj exp^iJ-nJ <S> l« aux ® C Wj8)Wi+1 (K) | 

7=1 
JV/2 



(16) 



containing only local operators for block Hilbert spaces. In equation ( jig ) we identified H-j-c I = 
&2i—i • <?2i + f?2i ■ (72i+i to separate the dependence on the parameter K. The block decomposition 
of the functional dependence O in ( jig ) contains two parts defined as 

eX H%W m ) I ( K )=<*Pl KH nA (17) 

and 

c, w J r ®??j+r<8.('Ha„x) f fJ _, ( K ) = c n!®H I+1 (K) (18) 



following the nomenclature for the block decomposition in the original worku. Ignoring the corre 



k@. 



lotion block part (|18| ) in (16) yields a decomposition of the quantum chain into commuting system 
block operators JlTj) analogous to the classical situation. This approximation is valid in the high 
temperature limit T — > oo where higher order terms of the coupling K, included in the correlation 
block part, vanish. _ . 

The correlation block part ([18]) is calculated using the Baker-Campbell-Hausdorf formulaE£l for 
blocks 

exp [K(H Hl + H Hl+1 )] = exp (K H Hl ) exp (K H Ul+1 ) ■ C Hl ®n I+1 (K) 
with C Hi ®h i+1 (K) = exp I — [H Hl , % +1 ] + •• • j • 

Relation © is an example for a pn^c* block decomposing. The separated correlation block part 
C?ii®Hi + i (K) of the functional dependence describes the quantum correlations between adjacent 
blocks in the decomposition. To eliminate the correlation block part ( |l8| ) in relation jlq ) an 
auxiliary space needs to be defined to describe the boundary conditions between adjacent blocks. 
Dependent on the choice of the auxiliary space (H a ux)/ the action of the block operators on the 
auxiliary space is determined. 

To describe the quantum correlations between the block I and the neighbouring blocks 1—1 and 
I + 1 the correlation block operator ( |l8| ) includes the coupling between the nearest neighbour 
(n.n.) single site spins of adjacent blocks. To include this n.n. coupling into the description of 
the boundary conditions we construct an auxiliary space including the n.n. spin sites of the block 
Hilbert space by choosing copies of the n.n. sites_pf the system block as visualized in figure pi 
According to the nomenclature in the original workd we call the constructed space shown in figure 



superblock 



<n i -2> 
aux 




n i +1 <n i +2> 
aux 



Figure 1. A system block enlarged by copies of the two nearest neighbour sites. The additional sites 
represent the auxiliary space and all single site Hilbert spaces together define the superblock. 

p] a superblock. To distinguish between the original system block sites and the copies of the n.n. 

sites in the superblock we mark the site indices of the copies with brackets '< . >'. 

Using the constructed auxiliary space together with relation (hq) results in the approximation 



Z 



xxx [On®Hw 
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{o\,, J odd} 



N/2 
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{/7-\systcm 
U Hi®(H a 



)/ K )} 



(19) 



Unlike in the classical case of section p] the block decomposition in relation (Il9|) does not allow 
for an exact conservation of the partition function. Furthermore, according to relation (|l9|), both 
auxiliary sites need to be truncated .within the GRRGT and by the choice of 7i aux we define an 
example for an active auxiliary spaceu. 

Here we give two remarks on the foregoing calculation: The choice of the particular auxiliary space 
allows for describing the boundary conditions between the blocks which determine the quantum 
correlations. The auxiliary space contains copies of the n.n. sites and neglects the effect of next 
nearest neighbour and further higher order couplings. Secondly we need to ensure that the auxiliary 
sites are treated as copies of the original sites during the GRRGT. Otherwise only the block Hilbert 
space Hi would have been enlarged and the description of the boundary conditions will fail. Since 
the single site Hilbert spaces and their copies are formally indistinguishable the identification as 
auxiliary sites in the superblock Hi (S> (H aU x)i is accomplished by the embedding and truncation 
operators G-h< $ mi ) and GT, „ , n s as illustrated in figure g. On the right hand side of figure 



superblock 



embedding 




Figure 2. The local embedding G-u'isn' an ^ truncation G^„ H procedure within the superblock. The 
lightly shaded dot is truncated within the spin decimation of the system block. 



the system block is visualized as the effective Hilbert space H' and the lightly shaded point 
denotes an even site which has been truncated in the GRRGT. 

After the necessary definitions within the GRRG approach have been performed the local RGT is 
given by the commuting diagramQ 



U\ 



'H> 



>n> 



i(*) 



H'j ^ 



% ®(«aux)l 



Hi ® (H aux ) 7 

Hi ® (W a ux)j 



(20) 



Due to the active auxiliary space no auxiliary sites are left after applying the GRRGT. Within a 
successive application of the GRRGT new copies of the changed n.n. sites for a system block need 
to be generated. 
We summarize the local operators as 

H/®(Waux) / ( K ) = CX P [^%8(W,ux)J 

= exp [K (<?(2i-2) ' &2i-l + <? 2l -l • <f-2i + &2i • CT2J+1 + <? 2l +l ■ <?(2j+2))] , 

Wl (K') = exp [K'H H , + f] 

= exp[K' (<J2i-i • (?2j+i) + /'] 



^H, ® (H,„), 



tr{cr 1: icvon} ® Ijcr,, iodd} ® tr H a ux 



with i G {!,..., iV/2} and 



with ? G {!,..., TV} . 



(21) 



By using (|l]) and calculating 

n tr <2*-2),(2i+2) exp [^i?H,«,(W aux )J 



°C^),( K )° 



G W/«>(«Lx)/ 



/ea 



tr e 



/ea 



= tr e 



> tr (cvcn) tr <cvon) II eXP [ K H Kii 



'(«««)/ J 
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(22) 

we proved that the GRRGT is perfectQ. In the final equation of d22) we have to trace over two 
copies of the even sites denoted as tr a and tr b . 

The choice of the auxiliary space does not allow for an exact treatment of the quantum correlations 
during the local RG procedure. Our approach therefore yields a perfect instead of an exact GR- 
RGT. Although the product decomposition ( J16| ) allows for further improvement in the description 
of quantum correlations by increasing the number of copied neighbouring sites it is not possible 
to construct an exact GRRGT. The definition of an exact GRRGT is possible but requires a dif- 
ferent and_more abstract auxiliary space. Results on the exact GRRGT will therefore be reported 
elsewherelia. 



IV. THE GRRG FLOW BEHAVIOUR IN THE FERROMAGNETIC REGIME 

To calculate the RG flow behaviour of the constructed GRRGT for the ferromagnetic isotropic 
spin-i Heisenberg chain we have to solve ( p0| ) for the effective coupling K' . It is convenient to 



rewrite the local operators in matrix form and solve the resulting set of equations. Using (10) the 
matrix representation of O-w (K 7 ) is given by 

Ou'j(K!) = a t 2i -i,2i+i + b' a 2l -i ■ &2i+i (23) 

and the coefficients are determined as 

a'(K',f) = [cosh^iO-sinh 3 ^')] e f ' and 

b'{K',f) = [sinh(^')cosh 2 (^:')-cosh(iY')sinh 2 (^')] e 7 ' . 

(24) 

Here we made use of the relation 

e K' (ff2i-i-a 2 i+i)+f _ e K' aZ t _ 1 <7% i+1 . & K' f|f 4 _ 1 «T* i+1 _ g A" v% i _ 1 <r% i+1 . g /' 

together with a trigonometric expansion using {<J2i-i a 2i+i) 2 — l2i— i,2i+i i a = x,y,z and a x ,a v 
and cr z as defined in (|l0|). According to ( |2l|) the functional dependence is invariant under an 
application of the GRRGT and we define 



G 



H, ® (Haux), °Hi®(H.^) I (K) = a(K) lai-l,2t+l + b(K) o 2 %-\ ' <? 2 j+l . (25) 



In the appendix we prove that this relation is well defined. Inserting (|23| ) and fl25|) into the local 
GRRGT (EOft and solving the resulting set of equations yields 

a(K) =a'{K',f) and b{K) = b' \K> ', f) . (26) 



To solve the coupled equations (|26|) the parameters a and b need to be calculated. Taking the trace 
over all the odd sites in ( f25| ) results in 

a(K) = lix{o Hl9lnm)l (K)}. (27) 



From (21) we conclude the diagonalizability of CVo®('H aux ) using a unitary transformation D = 
W HU with D a diagonal matrix. By identifying tr {O} — tr {exp [D] } we explicitly calculate 
b(K) as 

&(#) = ^ tr { 17+ ffai-i ■ $u+i 17 £/ f O U } 

= ^ tr { U ] CT 2 i-l ' <? 2l+ i t 7 ex P [- K '' D Wj®(Waux) J ] } ■ 

(28) 

Whereas the computation of the effective parameters a' and b' in (|2J) results from a reformulation 
of (Elf) the calculation of a and b involves the truncation procedure of the GRRGT. The matrix 



exp [ifO^gjH^) 1 used in the calculation of a and b is a diagonal matrix and the nonzero 
elements are the Boltzmann weights 



exp [KD Hli 



'(Wau*),] 



exp 



k^ Ej 



(29) 



of the superblock Hi <g> (W a ux)j where Ej denotes the corresponding energy eigenvalue. 
The core concepts of the GRRG follow the fundamental ideas of the density matrix renormalization 
group (DMRG) as explained in the original GRRG worku. Relation (|2Jj) defines the density matrix, 
originally introduced by R.P. FeynmanEj, of the superblock. The density matrix contains the 
information about the 'statistical importance' of each eigenstate of 0H /( g)(-H aux ) at temperature T. 
Compared to the numerical DMRG procedurcEH in which only the ground state (target state) 
is used for the calculation of the flow behaviour the GRRG uses all eigenstates weighted by their 
importance to define the local RGT. Within a finite temperature RG approach we expect that 
all the eigenstates of the superblock will contribute to the RG flow behaviour. The conventional 
strategy of using selected eigenvectors for constructing projection operators G-j-c and G^ „ , H ■. , 

used to project on a subspace of the total Hilbert spacecl, is therefore not sufficient in a finite 
temperature approach. 

For an explicit calculation of the RG flow, we have to solve ( J2q ) for the effective parameter K'(K). 
By means of (24) we obtain 



K'(K) = \ ■ In 
v ' 4 



a(K) + b(K) 



a(K) - 3b(K) 



(30) 



and a similar expression is obtained for the energy shift f'(K). Converting (p0[) into an explicit 
function of the flow parameter K requires the computation of the traces in (p7|)and ( |28| ) involving 
large matrix expressions. The resulting flow equation K'(K) displays a complicated structure and 
contains correction terms with increasing relevance in the low temperature regime. 
Figure shows a plot of the RG flow for the local RGT (EQ) using (El) in the ferromagnetic regime, 
i.e. K > and K' > 0. In figure g we also plotted the RG flow of the original system block with 
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Figure 3. The calculated GRRG flow for the superblock in the ferromagnetic regime (solid curve). The 
dashed curve displays a RG flow using fixed boundary conditions for the system blocks, only valid in the 
high temperature limit. 



no auxiliary space. This flow behaviour is equivalent to the one calculated by M. Suzuki et al 



§ 



valid in the high temperature limit. The different RG flows deviate from each other except in 
the high temperature limit, i.e. K — > 0, where the correlation block part (Iq) vanishes. From 
the curves shown in figure we conclude that the high temperature fixed point regime can be 
explored without defining an auxiliary space describing the quantum correlations within the block 
decomposition. However, figure [| also illustrates the necessity of including correlation block terms 
for computing the RG flow towards lower temperatures. 

According to the observed importance of the correlation block part by approaching lower temper- 
atures, we like to improve the local GRRGT including higher order correction terms. The special 
choice of the auxiliary space does not allow to describe quantum correlations beyond the nearest 
neighbour sites of the system block. However, enlarging the auxiliary space by including the next 
nearest neighbour (n.n.n.) sites results in the construction of an enlarged system block. Accord- 
ing to the local embedding and truncation maps (QW the additional single site spaces can not be 
marked as copies within the enlarged superblock. 

Instead of changing the embedding and truncation maps which in turn demands for choosing a 
different invariance relation ([L3|), we vary the size of the original system block, i.e. the scaling or 
reduction factor A in the GRRGT defined as 



A = 



number of lattice sites in the original spin chain 
number of lattice sites in the truncated spin chain 



(31) 



The previous calculations were based on a one site decimation procedure equivalent to a local 
GRRGT with a reduction factor A = 2. 

Figure displays two constructions of a superblock for a system block containing four single 
site Hilbert spaces with a reduction factor A = 3. In figure [| a superblock is defined with an 



enlarged superblock 
superblock 



Figure 4. A superblock and an enlarged superblock construction for a GRRGT with a reduction factor 
A = 3. The lightly shaded dots are decimated in the local GRRGT. 

auxiliary space including copies of the n.n. spin sites of the system block and an enlarged superblock 
containing also copies of the n.n.n. spin sites. According to the geometry of the system block the 
n.n. spin sites and the n.n.n. spin sites must be truncated within the local GRRGT which is 
consistent with our definition of the embedding and truncation maps for an active auxiliary space. 
The enlarged superblock contains four original system block sites and four further auxiliary sites. 
We define a goodness Q of the local GRRGT as the ratio of the number of copies of spin sites 
contained in the auxiliary space divided by the number of spin sites within the original system 
block 



Q = 



number of copies of spin sites in the auxiliary space 
number of spin sites in the system block 



(32) 



If no auxiliary space is defined Q = 0, whereas Q > 1 if the auxiliary space contains more copies of 
spin sites than original spin sites are contained in the system block. A sequence of improved local 
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GRRGTs is generated by enlarging the auxiliary space as visualized for a four site system block 
in figure 

In figure [| the RG flow for the superblock and the enlarged superblock structures depicted in 
figure are plotted. Again the dashed curve denotes the RG flow of the original four site block 
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Figure 5. The superblock RG flow and the enlarged superblock RG flow in the ferromagnetic regime. The 
dotted curve shows the flow behaviour using the enlarged superblock as compared to the usual superblock 
construction (solid curve). The dashed curve displays a RG flow using fixed boundary conditions for the 
system blocks, only valid in the high temperature limit. 

without any auxiliary space. All different RG flows plotted in figure ra show the correct high tem- 
perature flow behaviour by converging to the trivial high temperature fixed point. Apart from 
small corrections the superblock GRRGT and the enlarged superblock GRRGT display the same 
flow behaviour, indicating that the auxiliary space constructed by copies of n.n. sites provides 
a sufficient description of the quantum correlations in the plotted regime. However, we expect 
different flow behaviour for both superblock GRRGTs at lower temperatures. Figure ra shows the 
flow behaviour of the local GRRGT constructed from the superblock and the enlarged superblock 
displayed in figure away from the high temperature limit. 

Here we give a comment on the comparison between the GRRG method and the numerical DMRG 
procedure invented by S.R.WhiteEj which provided the core concepts for our workQ. The DMRG 
algorithm is designed for the numerical calculation of ground state properties of the physical sys- 
tem and is restricted to T = 0, i.e. no RG flow behaviour can be examineo£!rt3. Increasing the 
size of the system block in the DMRG method improves the numerical accuracy, but does not 
allow for studying the system at finite temperature. Within a DMRG calculation it is not possible 
to compute thermodynamic quantities and examine a phase transition at finite temperature. Ac- 
cording to the concepts of statistical physicsa a phase transition at finite temperature occurs as a 
nontrivial fixed point in the RG flow whereas the high and low temperature limits represent trivial 
fixed points regimes. The RG flow behaviour in a nontrivial fixed point regime is characterized by 
a set of critical exponents defining the type of the phase transitioned. 

From the GRRG method we calculate the thermal critical exponent v representing the magnetic 
phase transition usingca 



k)g( 



dK' 
8K 



K< 



log(A) 



(33) 



where K* denotes the value of the fixed point coupling in the RG flow equation (BOf). ._. 
Strongly correlated systems do not exhibit a nontrivial fixed point in one dimension!!^!. Solving 
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Figure 6. The RG flow of the enlarged superblock GRRGT (dotted curve) in a lower temperature regime 
compared to the GRRG flow behaviour of the superblock (solid curve). 



the RG flow equations plotted in the foregoing figures therefore leads to the same trivial flow 
behaviour as in the classical analogue examined in section [En without exhibiting a nontrivial fixed 
point. Nontrivial fixed points occur in quantum spin chains of higher dimensions. 
In 1975 A. A. Migdal proposed a method for analytical continuation to highes, dimensions of RG 
recursion formulas for strong coupled systems exhibiting global symmetriesc§c£L The result of 
A. A. Migdal, applicable to a variety of decimation and truncation proceduresEO, was rederived 
and rigorously analyzed by L.P. Kadanoff by inventing the Kadanoff bond moving proceduraa. 
Both authors assumed a model Hamiltonian with n.n. interactions. Although L.P. Kadanoff has 
generalized and corrected the results of A. A. Migdal, the resulting formula for isotropic quantum 
spin models was exactly the same as the equation proposed by A. A. Migdal given by 



K(XL) = X d - X R x [K(L)} 



(34) 



In relation ( p4j ) L denotes the lattice constant, i.e. the distance between two n.n. spin sites and 
the functional R x denotes the RGT in the coupling K. According to the notation used by L.P. 
Kadanoff XL denotes the lattice constant in the decimated spin chain. The calculations presented 
in this work include no explicit dependence on a lattice constant L and we identify K' = K(XL). 
We applied the Kadanoff bond moving procedure to the isotropic XXX spin-i Heisenberg model 
in dimension d = 3 which exhibits a nontrivial fixed pointO. We determined the nontrivial fixed 
point for all constructed local GRRGTs and confirmed that in dimension d < 3 all GRRG flows 
exhibit only the trivial fixed points K* = and K* — oo. 

In table | we have summarized computed fixed point values together with the corresponding crit- 
ical exponents v calculated by equation (f33|) for dimension d = 3. We applied the GRRGT with 
a reduction factor A = 3 and A = 4 using the superblock and the enlarged superblock struc- 
ture. The calculated critical exponents vary between 0.47 and 0.49 ordered by the goodness Q 
of the approximation used for describing the boundary conditions. For both reduction factors we 
furthermore calculated the critical exponents if no auxiliary space is provided corresponding to 
an approximation only valid in the high temperature limit. The calculated values of the critical 
exponent deviate significantly from the results obtained by using the superblock or the enlarged 
superblock structure. In table 1 we further compared the results of the GRRG method with the 
outcome of other methods from the literature. The values of the critical exponent for the Approx- 
imate decimation method and the Mean Field RG method deviate significantly from each other. 
Combining both methods resulted in an even higher value for the critical exponent as compared 
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Method 


RG flow fixed point K* 


Critical exponent v 


Goodness Q 


A = 3 reduction and no auxiliary space 


0.522 (±1.0- KT 4 ) 


0.645 (±1.0- 10" 4 ) 





A = 3 reduction superblock 


0.640 (±1.0- 10^ 4 ) 


0.470 (±1.0- 10~ 4 ) 


0.5 


A = 3 reduction enlarged superblock 


0.627 (±1.0- 10~ 4 ) 


0.489 (±1.0- 10~ 4 ) 


1.0 


A = 4 reduction and no auxiliary space 


0.703 (±1.0- 10~ 4 ) 


0.595 (±1.0- 10" 4 ) 





A = 4 reduction superblock 


0.851 (±1.0- 10~ 4 ) 


0.469 (±1.0- 10" 4 ) 


0.4 


A = 4 reduction enlarged superblock 
Approximate decimation methodli3 
MFRG combined with decimationu 
Mean Field RG (MFRG)B 


0.837 (±1.0- 10~ 4 ) 


0.475 (±1.0- 10~ 4 ) 


0.8 


0.344 


0.714 





0.312 


0.758 


- 


0.275 


0.450 


- 



Table I. The numerical fixed point values and the corresponding critical exponents u of the isotropic 
quantum spin-^ Heisenberg model in dimension d — 3 calculated by different methods. 



to the approximate decimation method, difficult to validate. For strong coupled quantum spin 
lattices at finite temperature T no exact approach is available in dimension d = 3 to calculate 
thermodynamic behaviour. The GRRG method with the proposed auxiliary space is a rigorous 
and analytic approximation. The approximation is controlled by a goodness parameter calculated 
in equation (133) yielding consistent results. 



V. THE ANTI-FERROMAGNETIC ISOTROPIC HEISENBERG CHAIN 

In this section we examine the anti-ferromagnetic regime, i.e. K < 0. By using a reduction 
factor A = 2 or A = 4 in the local GRRGT the anti-ferromagnetic part of the RG flow exhibits an 
unphysical behaviour, i. e. applying the local GRRGT once yields a ferromagnetic coupling K > 0. 
The situation is different for a system block structure with a reduction factor A = 3. According 
to the geometry of the enlarged superblock the GRRG flow shows the correct anti-ferromagnetic 
behaviour. 

Due to the inherent global symmetry of the isotropic quantum spin-^ Heisenberg model the eigen- 
states of the model Hamiltonian can be represented by the spin z-component for each lattice spin 
site. Using this representation the ground-state for the quantum spin-^ Heisenberg model is rep- 
resented by an alternating sequence of spin up and spin down z-components. Figure |7] visualizes a 
RG-step, i.e. applying the RGT once, using a reduction factor A = 2. The original n.n. coupling 
of anti-ferromagnetic type is transformed into a n.n. coupling of ferromagnetic type according to 
the structure of the system block. By using a reduction factor A = 2 the geometry of the spin 
lattice does not allow for calculating an anti-ferromagnetic GRRG flow behaviour. 
By choosing a reduction factor A = 3 the geometry of the system block changes. In figure @ we 
depict a RG-step for a system block composed of four sites. By applying the RGT the original 
anti-ferromagnetic n.n. coupling K < is transformed into an effective anti-ferromagnetic n.n. 
coupling K' < 0. According to the structure of the system block it is possible to construct an 
entire anti-ferromagnetic GRRG flow using the enlarged superblock structure. 
In figure the anti-ferromagnetic GRRG flow using the enlarged superblock structure depicted 
in figure |g is plotted by the dotted curve. The dashed curve displays the flow behaviour using a 
system block containing four spin sites without an auxiliary space. Although the RGT respects the 
geometry of an anti-ferromagnetic quantum spin chain no quantum block correlations are taken 
into account which is valid in the high temperature limit T — > oo. In figure M two further GRRG 
flows were plotted using a superblock structure with a reduction factor A = 2 and A = 3. The 
lighted solid curve shows the expected unphysical behaviour for the reduction factor A = 2 as de- 
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Figure 7. An application of the RGT for a reduction factor A = 2. The nearest neighbour (n.n.) coupling 
K < is transformed into an effective n.n. coupling of ferromagnetic type K' > 0. 



picted in figure J7|. The black solid curve displays the GRRG flow according to the superblock with 
a reduction factor A = 3 depicted in figure 0. Although the system block, containing four original 
spin sites, respects the anti-ferromagnetic geometry, the auxiliary space does not. However, in the 
high temperature regime the GRRG flow for the superblock and the enlarged superblock structure 
using a reduction factor A = 3 display equal behaviour. In the finite temperature regime ground 
state properties become increasingly important and the auxiliary space for the superblock with a 
reduction factor A = 3 does not provide the correct anti- ferromagnetic boundary conditions. 



VI. CONCLUSIONS 



In this paper we applied the Generalized Real-space Renormalization Group (GRRG) method 
to the isotropic spin-^ Heisenberg model. The analytic method is non-perturbative and yields an 
approximate RG flow behaviour in the finite temperature regime controlled by a quality measure, 
i.e. the goodness parameter. 

In section [II we constructed a local GRRG transformation based on a chosen auxiliary space. The 
spin site decimation was formulated in terms of an embedding and truncation map. 
In section IV we examined the ferromagnetic RG flow behaviour for the superblock and the enlarged 
superblock structure using different reduction factors for the spin site decimation. We explored 
the RG flow behaviour in three dimensions using the Kadanoff bond moving procedure originally 
introduced by A. A. Migdal as a method for analytical continuation to higher dimensions in RG 
recursion formulas. The nontrivial finite temperature fixed point of the different RG flows was cal- 
culated and the related critical exponent was examined according to the goodness of the auxiliary 
space approximation. The critical exponents computed by the GRRG method were compared to 
other results from the literature, mostly calculated by numerical techniques without providing a 
quality measure of the approximation. 

In section |V| the anti-ferromagnetic part of the flow behaviour was explored. Only system block 
structures with an odd reduction factor allowed for a local GRRG transformation respecting the 
anti-ferromagnetic geometry of the quantum spin chain. The enlarged superblock GRRG flow for 
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Figure 8. An application of the RGT for a reduction factor A = 3. The anti-ferromagnetic nearest 
neighbour (n.n.) coupling is transformed into an effective n.n. coupling of anti-ferromagnetic type K < 0. 

a two site decimation in the quantum spin chain exhibited the correct physical flow behaviour 
and was proved as a valid approximation in the finite temperature regime. The auxiliary space 
of the superblock for the two site decimation did not provide the correct description of the anti- 
ferromagnetic boundary conditions resulted in an unphysical flow behaviour in the lower temper- 
ature regime. 

Applications of the GRRG method to other quantum systems such as the Hubbard model will be 
reported in the future, although further development of auxiliary spaces is the core concept of the 
GRRG approach. In this work we showed that away from the high temperature fixed point regime 
boundary conditions become increasingly important in the calculation of the RG flow behaviour. 
An exact GRRG method provides all possible boundary conditions at an arbitrary temperature. 
We therefore have constructed an exact local GRRG transformation for the spin-i Heisenberg 
model using a passive auxiliary space providing all necessary boundary conditions between adja- 
cent system blockala. 

Although the GRRG method is designed as an analytic RG approach we are working on a numer- 
ical implementation of a GRRG algorithm applicable to nonlinear partial differential equations. 
A detailed description of the mathematical formulation and the numerical implementation of the 
algorithm will be presented in the near future. We compute critical exponents to determine the 
universality class of the nonlinear partial differential equation with a reduced number of degrees 
of freedom. 
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Figure 9. The lighted and the black solid curve display the unphysical GRRG flow behaviour using a 
superblock structure with a reduction factor A — 2 and A = 3 respectively. The dotted curve depicts the 
anti-ferromagnetic GRRG flow behaviour for the enlarged superblock with a reduction factor A = 3. The 
dashed curve displays the flow behaviour for the system block with a reduction factor A = 3 using no 
auxiliary space. 

Appendix 

In order to derive equation ( pq ) we introduce the abbreviations 

tr OV cn {exp [KH ni ^ naux)i ] } = tr cvcn { ¥?({<?*•},£/) } = x(?2i-i,?2i+i) 



and <p and \ are functions of the Pauli spin matrices defined in ( |10[ ) . Using the rotation symmetry 
we derive an equivalent representation {(?;} by 



?{&}) = ?({%}) 



with the rotation map R TOt defined as 



From 



we deduce 



X(<?2i-i,<?2i+i) = tr even | <p({(?i}) | 
= tr cvcn / J <p({ai}) | 
= tr even / J ip({a' t }) J 

= x(4i-i)4i+i) 

X(a 2 i-l,(T2i + l) = X(^2i-l ' CT 2l+ l) . 
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Using (p3|) we derive (p5|). 



(i 
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